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The application of matrix iterative analysis to the solution oj waveguide 
discontinuity problems is discussed. It is concluded that the "Gauss-Seidel" 
or "point-single-step" method offers several advantages over more conven- 
tional invertive procedures, particularly in the speed oj execution. Two 
examples are presented as illustrations: analysis oj an H-plane discon- 
tinuity in a rectangidar waveguide and conversion jrom TE U to TM U 
modes at an abrupt discontinuity in a circular waveguide. The latter 
results are shown to be in good agreement with measured values obtained 
in a previous investigation. 

I. INTRODUCTION 

The analysis of waveguide discontinuities, for application to the 
design of antennas and microwave networks, continues to offer challeng- 
ing problems in electromagnetic theory and microwave engineering. 
Thus far, the solution of these problems has depended to a large extent 
on various approximate techniques, such as variational and quasi- 
static methods, 1 which are extremely useful but nevertheless limited 
in applicability. 

The shortcomings of classical analysis have been surmounted to a 
large extent by our ability to solve electromagnetic boundary value 
problems by numerical methods, making extensive use of digital com- 
puters. Computational techniques are not only an abundant source 
of engineering data, which might otherwise require elaborate construc- 
tion and experiment, but they can also provide a unique analytical 
laboratory in which to evaluate approximate theoretical methods under 
easily controlled conditions. In this paper, we shall be concerned with 
these numerical methods as they apply to certain waveguide discon- 
tinuity problems. 
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For the sake of simplicity we shall consider, as an example, the 
problem of two waveguides with similar cross sections connected to- 
gether at the plane z = 0, as illustrated in Fig. 1. A wave is shown 
incident from the smaller waveguide impinging on the discontinuity. 
The result will, of course, be to excite an infinite number of normal 
modes in each guide, some of which carry real power away from the 
junction, with the remainder being evanescent and contributing to 
the electromagnetic field only in the vicinity of the connecting aperture. 
It must be recognized that these evanescent modes play an important 
role since they, in part, determine the amplitudes and phases of the 
propagating modes. It is the fact that an infinite number of waves 
must, in principle, be considered that makes this type of problem so 
difficult. 

The contents of the paper may be summarized as follows: We begin 
by establishing an appropriate form of the uniqueness theorem for 
Maxwell's equations as they apply to boundary value problems of 
this type. In numerical analysis, the criteria for uniqueness are of more 
than academic interest since they provide meaningful and practical 
methods by which to assess the accuracy of results. Next, the normal 
mode representation of the fields is discussed, the object being to 
arrive at a matrix equation formulation of the problem in which the 
components of the unknown vector are the modal coefficients. It is 
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Fig. 1 — Waveguides of similar cross section connected at the plane 2 = by 
an abrupt discontinuity. A wave is assumed incident from the smaller guide. 
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suggested that this matrix equation may be solved by an iterative 
procedure and, upon studying the convergence properties of such 
methods we find a critical dependence on the particular algorithm 
used. Two examples will be presented as illustrations, analysis of an 
H-plane discontinuity in a rectangular waveguide, and conversion from 
TE n to TMn modes at an abrupt discontinuity in a circular waveguide. 
The latter results are shown to be in good agreement with measured 
values obtained in a previous investigation. 

Rationalized MKS units and the (suppressed) harmonic time de- 
pendence exp (—iut) will be used, unless otherwise specified. 



II. UNIQUENESS AND ERROR CRITERIA 

A representation of the discontinuity is shown in Fig. 2. It is assumed 
that the regions to the left (denoted by — ) and to the right (denoted 
by +) are each filled with homogeneous material, but with possibly 
different constitutive parameters. Maxwell's curl equations in the 
respective regions are thus given by 



V X E* = -V 



air 

dt 



(1) 



V X IT = e a 



dt 



As usual for uniqueness theorems, we begin with two solutions in 
each region presumed to be correct, and denote the differences respec- 
tively by E*, H*.* Then from the Poynting theorem, 2 it follows that 
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Fig. 2 — Waveguide discontinuity showing boundary surfaces A, B, C, D and 
respective normals ni, n 2) n 3 , n.|. 



* Physically, these fields would correspond to a waveguide discontinuity problem 
without excitation. 
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in the — region 

jf (E7 X H~) n, dS + J J (E7 X HI) n 2 dS 

- -" s //£ E ~ E ~ dV - "" £ //t H ~ H ~ rfF (2) 

and in the + region that 
jj (El X H+)-n 3 d£ + // (El X H>n 4 tfS 



6 

5/ 



where d/dt denotes differentiation with respect to time and the sub- 
script t denotes the field transverse to the generatrix of the cylinder. 
The unit normal vectors ni , n 2 , n 3 , and n 4 are shown in Fig. 2. 

One must also take into account the fact that certain physical con- 
siderations will limit the class of admissible solutions. For example, 
if we let the surfaces A and D recede to infinity, then all evanescent 
modes will have decayed to zero and the respective surface integrals 
then represent the power flow away from the discontinuity. Assuming 
no loss, the total power must vanish. Furthermore, it can be shown 
from Maxwell's equations that the transverse components of electric 
and magnetic field at the interface must be continuous. Adding (2) 
and (3), we find that the following time derivative must vanish, 



at |_ € JJJv- 



E--E - dV + yT fff H~-H~d7 

+ 6 + fjf + E + E + dV + m + /// H + -H + dVJ = 0, (4) 

We may, however, regard the quantity in brackets as having had 
a zero value at some time, say at t = 0, the excitation time. The term 
in brackets therefore, vanishes for all time and, since each of the in- 
tegrands is positive semi-definite, they must vanish separately. Thus, 
at each point in the + and — regions, 

Ej — E 2 = Ha — H 2 = ,-x 

E7 - E7 = Hr - Ho = 

and the solution is thereby shown to be unique. We may now state 
the following uniqueness theorem for waveguide discontinuity problems. 
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Theorem: The solution to a waveguide discontinuity -problem is uniquely 
specified if it can be shown to have the following properties: 

(j) It satisfies Maxwell's equations and the appropriate boundary 
conditions in the regions on each side of the discontinuity. 

(ii) The components of electric and magnetic fields tangent to the 
interface are continuous. 

(Hi) In the case of a lossless discontinuity, energy is conserved. 

These three conditions obviously play an important theoretical role 
in the solution; where numerical methods are used, they also provide 
fundamental criteria by which the accuracy of computed results can 
be assessed. Accordingly, we shall define the following quantities to 
be used as error criteria: First, there is the parameter e p , which indicates 
how well the solution conserves energy, given by 

Cp = —B 1. (6) 

where P T ,P t , and P iD0 are the reflected, transmitted, and incident powers, 
respectively. Second, the mean square error in the tangential electric 
field is defined by 

[[ | (E; - E~) | 2 dA 

e s = JJApert,ir " (7) 

[f | E! ino) | 2 dA 

J "Aperture 

and third, for the magnetic field, 



!L 



I (a: - K~ t ) | 2 dA 

Cu = ^^ J (8) 

«*Apertura 

where E Cino> and H (,no> refer to the incident wave. 

The smaller the quantities e p , e E , and e h , the more closely the 
boundary conditions are satisfied, at least in the mean square sense, 
and the more accurate we shall consider the solution to be. 

III. MATRIX FORMULATION OF THE BOUNDARY VALUE PROBLEM 

The most convenient format for numerical solution of waveguide 
discontinuity problems is a matrix representation, in which the modal 
coefficients form the unknown column vectors and the discontinuity 



654 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1967 

is characterized by a square matrix. We recognize that there is also 
an analogous integral equation in terms of the aperture electric or 
magnetic field. However, since the numerical solution of the integral 
equation is generally carried out by reducing it to a matrix equation, 
we shall proceed to the matrix formulation directly from the physical 
characterization of the boundary value problem. This matrix equation 
will then be solved by an iterative method, the theory of which is 
discussed in Section IV. 

It is assumed that in each of the waveguides, the electromagnetic 
fields may be characterized by a denumerable set of known vector 
eigenfunctions which may be ordered according to some index. We shall 
be concerned only with the transverse fields,* denoted as follows: 

+ E£(r) (p => 1, 2, 3, • • •) denotes the transverse electric field for the 
pth TM mode in the + waveguide, with r as the position vector in 
the transverse plane. 

+ H£(r) = transverse magnetic field for the pth TM mode in the 
+ waveguide. 

+ E£'(r) = transverse electric field for the pth TE mode in the + 
waveguide. 

+ H£'(r) = transverse magnetic field for the pth TE mode in the 
+ waveguide. 

By replacing the + by — we have the analogous notation for the 
other waveguide. An important point concerning sign convention is 
that the unknown modes in the — waveguide will all be taken to 
propagate away from the discontinuity, i.e., in the —z direction. 
Although the electric field does not change sign when the direction 
of propagation is reversed, the magnetic field does, and this fact must 
be carefully taken into account. 

In order to define the amplitudes of the respective vector wave 
functions, we adopt the following normalization/ written in terms of 
integrals over the waveguide cross sections: 

[ *E£-*E:* dA = |*A$| 8 *„ (9) 

EL'-*B"*dA = «V 5 P< , (10) 



./: 



in which h p is the respective characteristic wavenumber, and n is the 
permeability, which in our case will be the permeability of vacuum, 

* It is assumed that the individual waveguides can support pure TE and TM 
modes, which is the case for applications of interest here. 
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since both waveguides will be assumed empty. t By introducing the 
Kronecker delta 8 pq , we have also expressed the fact that the trans- 
verse fields in the individual waveguides are orthogonal. 

Once the normalizations for the electric wave functions are defined, 
those for the magnetic field are also specified since, for both the TE 
and TM modes, the transverse electric and magnetic fields are uniquely 
related. In particular, for a TM mode 

*h; = ±~ e, x *e; (ii) 



and for a TE mode 



*H£' = ±^e s X *E£' f (12) 



CO/J 



where e 2 is a unit vector in the z direction. Note again that the sign 
convention is such that a field in the — waveguide is taken to be a 
reflected wave, travelling away from the discontinuity. The magnetic 
field normalization is thus given by 



./: 



-R' p -*W*dA =«V« pt (13) 



± -a v ,± -a f q '* dA = | *#' | 2 5 P9 . (14) 

'±4 



Both sets of transverse wave functions have the property of com- 
pleteness, which is to say that any transverse electric (or magnetic) 
field can be synthesized from a set of TE and TM vector wave functions, 
provided that the directions of propagation of the normal modes are 
known. For the problems to be considered here, this latter information 
is available from physical considerations, since all modes propagate 
away from the junction with the exception of the incident wave whose 
amplitude is known. This amplitude will be taken to be that of a 
normalized mode. 

We now derive the appropriate matrix representation for the dis- 
continuity problem. Assume a dominant (TE) mode wave (E", H") 
is incident from the — guide, setting up a transverse electric field in 
the aperture just to the left of the junction. This field, referred to as ~E, , 
may be synthesized as follows : 

"E, = "EC + £ -A p -E' v + Z -B Q ~E' Q ' (15a) 



t The asterisk (*) denotes the complex conjugate. 
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with the modal coefficients ~A P and ~B Q as yet undetermined. The 
corresponding transverse electric field on the -f- side, denoted by + E, , 
would then be given in terms of normal modes on the + side by 



+ e ( = z + a p + e p + 2 + VEr. 



(15b) 



Since the transverse electric field is continuous across the aperture, 
we have that 

E ( = ~E, on C 
+ E« = on D-C. 

As shown in Fig. 2, D-C represents the conducting wall which makes 
up the remainder of the junction, and on which the transverse electric 
field must vanish. Expanding (15a) in a Fourier series of modes in 
the + waveguide, we find that the modal coefficients are related by 

+ A P = j^rrs ( -E('- + E£* dA + t-^tts Z ~A q f ~E' a - + E p * dA 

\ «j» | Jc ! "j> | <j = l Jc 

+ mTt-,2 Z ~B Q f "EJ'. + Ei* dA (16a) 
*B v = -h f -E['- + E p '* dA + -^it~A Q f ~E' a - + E p '* dA 

CO H Jc CO fi , = i Jc 

+ -h £ ~B m f -E' Q '- + E' P '* dA (16b) 

co n e= , J c 

or, more succinctly, using partitioned matrix representations, 



r a 



+ £>' 
-5-5 tf 



L^aJ 



+ 



2 2 «'J 2 2 ^ 

|_OI /Li CO /X 



[s r 


( — 
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8 r 


[— 


+ 
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i 




// 


i 




[~a~] 


8 r 


[— 


+1 


e r 
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+~ 




_~(B_ 
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// 




// 


n 







(17) 
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in which the vectors and submatrices are defined as follows: 




(IF,), = j ~E' t ,+ K*dA 



(18) 



(19) 



*SD' = 



1 


|2 





••• 


r« 







1 


••• 


1 + AJ I s 













(20) 



and d is the identity matrix. The matrix £, whose transpose appears 
in (17), is the matrix of coupling coefficients, defined as the scalar 
products of electric transverse vector wave functions for the wave- 
guides on each side of the discontinuity. The four index notation is 
interpreted as: 



-, + 



-/. 



E' • + E"* dA 



(21) 



with analogous definitions for other combinations. 

The system of equations given in (17) is clear[ 
since the number of unknowns is twice the numbed 
ever, an additional set can be derived by empl^ 
condition that the transverse magnetic field must 
across the interface. The matrix equation, ant 
corresponding to this second boundary condition, ij 



underdetermined 
)f equations. How- 
ling the boundary 

jo be continuous 
50US to (17), but 
pven by 



r ~a 




r°i 
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4m $ 

to E 


_~«j 




-9J 




_~SD" 
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2 2 

to e 



©' 



-J" 



3e 3 



JC 



+, - 



+, - 



oc' 



(B + 



(22) 
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in which 



9 = 



(23) 



2D" = 



"A." 




(24) 



and the matrix 3C, whose transpose appears in (22), is the matrix of 
scalar products of magnetic transverse vector wave functions. The 
four-index notation is interpreted in the same way as in (21). 

It should be noted once again that all the matrices which appear 
in (17) and (22) are infinite matrices, corresponding to the fact that 
in general an infinite number of modes are excited in the neighborhood 
of the discontinuity. In practice, of course, there must be a truncation 
and the problem then becomes one of solving a set of matrix equations 
whose order, N, ^depends on the accuracy required. Unfortunately 
there is, as yet, no iray in which the number of modes required to produce 
a given accuracy can be predicted. We can only emphasize the need 
for meaningful er»r criteria which will act as a guide in choosing a 
number of modes much will be large enough to give sufficiently accurate 
results but at the^une time not be so large as to require excess com- 
putation. It is exmbted that the criteria given in Section II will prove 
very useful in thiflwespect. 



IV. MATRIX ITI 

It was shown 
problem of int 
linear algebraic 
physics, so that 
of matrix equatioJ 



METHODS 

|e previous section that the waveguide discontinuity 
lere can be formulated in terms of a system of 
tions. This is a recurrent theme in mathematical 
tensive theory concerned with the efficient solution 
las evolved. In this section, we shall be concerned 
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with some of the elements of this theory, placing particular emphasis 
on the solution of matrix equations by iteration. 

The system of linear algebraic equations which results from satisfying 
the aperture boundary conditions on the transverse electric and mag- 
netic fields can be written in the matrix forms 

+ a = 01 + <R'~a (25) 

~a = V + S- + a, (26) 



where 






(27) 



the vectors 01 and V and the matrices (R and S being correspondingly 
identified from (17) and (22). Equations (25) and (26) are easily un- 
coupled to give 

[if - (RS]- + a = 01 + (TO (28) 

[si - S(R]-"a = V -f- SOL (29) 
both of which are seen to have the general form 

311a: = y. (30) 

In (30), x is an iV-dimensional complex vector whose components 
are the coefficients of the normal modes in the two waveguides, 3TI is 
an N X N complex matrix characterizing the discontinuity, and y is an 
excitation vector due to the incident wave. 

The obvious method of solving (30) is to compute the inverse of 3H 
and thus directly obtain 

x = an - y (3i) 

However, we should recognize that it is the solution vector x which 
is required, and that computing the inverse is not always the best equa- 
tion solving technique. For example, because the modal coefficients 
may decrease slowly with mode index, an accurate approximation of 
the physical problem often requires that 911 be a very large matrix, 
and inversion procedures for large complex matrices require considerable 
computational effort. An alternate approach is therefore suggested, 
namely the solution of (30) by a method of iteration. 

In an iterative algorithm, we begin with an initial "guess" for the 
solution and, from this, generate a supposedly improved solution, 
repeating the process until successive iterations give results which 
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-TE I0 MODE INCIDENT 

Fig. 3 — H-plane discontinuity in a rectangular waveguide. The incident mode 
is a TEio mode (electric field vertical). 

agree to within some prescribed norm. The solution to which the 
procedure converges must, of course, be independent of the initial as- 
sumption. 

A tempting iterative procedure for the present problem is suggested 
by writing (28) in the form 



+ a = 11 + <RT) + (RcS- + a 
with an initial assumption 

V" = ii + aru. 



(32) 



(33) 



Physically this corresponds to first assuming the aperture electric 
field to be that of the unperturbed incident wave, and calculating the 
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Fig. 4 — Transmission coefficient of an H-plane discontinuity. 



WAVEGUIDE DISCONTINUITY 



661 



TE and TM modal coefficients in the + waveguide on this basis. The 
corresponding magnetic field is then determined on the + side of the 
aperture and, with the aid of the appropriate continuity condition, 
is used to find the magnetic field and subsequently an "improved" 
electric field in the — region. This second guess for the aperture 
electric field is then used to repeat the process, etc. 

As a test, this algorithm was applied to analysis of an H-plane dis- 
continuity in a rectangular waveguide, the incident wave being a 
TE 10 mode of normalized amplitude [see (14)] as in Fig. 3. The dimen- 
sions were ka 2 = 4.5 and fca, = 3.5 where k is the free space wave- 
number. With this choice of parameters, only the TE 10 modes can 
propagate in each guide. (This problem is discussed in further detail 
in Section V.) 

Fig. 4 shows the result of calculating the real part of the modal 
coefficient for the TE 10 mode in the larger waveguide, plotted as a 
distribution of points giving the value at each iteration. The Fourier 
series for this particular example was truncated after twenty five terms. 
Aside from a small amplitude oscillation of less than one percent rms, 
the results seem reasonable, especially in view of calculations for the 
mean-square errors e E and e„ , which are illustrated in Fig. 5. These 
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Fig. 5 — Mean square errors in transverse electric and magnetic fields for an 
H-plane discontinuity. 
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Fig. 6 — Oscillatory instability of a nonconvergent algorithm for the H-plane 
discontinuity analysis. 

decrease monotonically with succeeding iterations, e b approaching 
approximately 0.001 and e h a value of about 0.012. The larger error 
in the magnetic field can be attributed to the singular behavior in 
H near the corner of the discontinuity. The asymptotic value of the 
energy parameter e P is approximately 0.007. 

The apparently accurate results obtained using this algorithm are 
in fact quite deceptive, and may actually be attributed to the propitious 
initial choice for the aperture electric field. It will be recalled that a 
very important criterion for validity of an iterative procedure is that the 
results be independent of the initial assumption. In order to determine 
whether such a criterion is satisfied for this particular algorithm, the 
TE 10 modal coefficient for the larger waveguide was arbitrarily doubled 
after the fourth iteration, which is equivalent to deliberately assuming a 
poor initial choice for the aperture field. The effect, shown in Fig. 6, 
indicates that the algorithm does not relax to the previous values, but 
continues to oscillate with a large amplitude. Similarly large fluctua- 
tions occur in e E , e H , and e P , the conclusion being that this particular 
precedure is not satisfactory, and can be expected to give reasonable 
results only if the initial choice is a very good one. The reason for this 
instability will become apparent after we consider those aspects of 
matrix-iterative analysis which are relevant to these problems. 

In the usual framework for iterative procedures, the matrix equation 
satisfied by the unknown vector x is written in the form 

x = mx + /, (34) 

where 3TZ and / are appropriate to the particular scheme being used. 
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This leads very naturally to the recursive formula relating the m + 1 
to the ra iteration, 

x lm+1) = mx (m) + /. (35) 

We denote the error vector at any iteration by £ <m) , where 



£ o») = ,.(«) 



- x, (36) 



and x is the exact solution to (34). Then, by substituting (36) into (35), 
we find that £ <m+1> is related to e im) by 

f(m+ „ = gil£ < m ) j (37) 

Therefore, the error at the rath iteration is expressible in terms of the 
initial error by 



£ 



(•"> ^ cn7 m „(°' 



= 9lTV 0) . (37a) 



For an absolutely converging solution we thus require that 

£ (m) — > as m — » oo (38) 

regardless of the initial guess x m . This is equivalent to 

3lT -» as m -> oo (39) 

where the in (38) and (39) denotes a null vector or matrix, respectively. 
It can be shown 4 that an N X N complex matrix 311 is "convergent", in 
the sense of (39), if and only if all the eigenvalues X,- of 971 magnitude 
less than unity, i.e., 

| X,- | < 1 all i. (40) 

We can easily see why this requirement will guarantee convergence, 
at least for the special case where the eigenvectors a,- of 9TC span the 
space of N-dimensional complex vectors. The initial error is then 
expressible as 

e w = Z On , (41) 

where the C,- are constants. The error at the rath iteration then be- 
comes, from (37), 

£<-> = f; c.-snTa,- = E CXbn (42) 

and, as m — > » f each term in the sum approaches zero, provided, 
of course, that (40) is satisfied. 
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Equation (42) also reveals useful information concerning the speed 
of convergence, which is seen to depend on how close the magnitudes 
of the X,- are to unity. Clearly if the largest eigenvalue is very near 
one in magnitude, the convergence will be very slow and hence a large 
number of iterations will be required. It would be logical, in the light 
of this reasoning, to evaluate the magnitude of the largest eigenvalue 
for the H-plane discontinuity problem discussed proviously. Unfor- 
tunately, because of the geometrical asymmetry, the matrix 3H is not 
Hermitian and so the usual computational techniques for determining 
eigenvalues cannot be used. We can, however, find an upper bound 
for the modulus of the maximum eigenvalue, p(9TC) = max {| X,- |}, 
given by 5 

p(l) ^ [p(9n: + 3TC:)] 1 . (43) 

where f denotes the conjugate transpose matrix. Note that 311*311 is 
Hermitian, so that standard computer programs can be used to evaluate 
its eigenvalues. We find for the previous H-plane problem that | X,- 1 ^ 1.16 
which, although, not conclusive, shows the possibility of such an os- 
cillatory instability. 

One technique which is suggested as a means of obtaining a con- 
vergent algorithm is called the "Richardson" or "point-Jacobi" method. 6 
In this approach, the matrix 91Z is first partitioned as 

911 = £> + £ + 11, (44) 

where SD is a diagonal matrix containing the diagonal terms of 311, £ 
is a strictly lower triangular matrix and 11 is a strictly upper triangular 
matrix. The system (34) is then written as 

(d - £>).r = (£ + 11)2 + / (45) 

from which 

x = (d - £>)"'(£ + m).T + (8 - £)) -1 / 

= m R x + /„ , (46) 

(46) being the matrix representation of the "Richardson" or "point- 
Jacobi" method. 

A modification of this procedure is referred to as the "Gauss-Seidel" 
or "point-single-step" iteration method. 6 Note that in solving (46) by 
iteration, the components of z lm+1) are all computed from the com- 
ponents of x lm) . Intuitively, it would seem more attractive to use the 
latest estimates of x, i.e., in computing x] m+l) we should use, wherever 
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they appear, the components xJT* 1 * (k < j) already computed, and 
in this way utilize the most accurate information available. This pro- 
cedure is, in fact, easier to implement in a computer program and, 
in addition to requiring less storage, often has better convergence 
properties than Richardson's method. It may be shown that the matrix 
representation, analogous to (46), for the "Gauss-Seidel" method is 

x = mi x + f G (47) 

3R = (tf - 2D - cCr'll (4g) 

fa -(*-»- £)-'/■ 

3), £, and 11 having been defined previously. 

The eigenvalue condition given in (40) is, of course, a very restrictive 
one, so that iteration procedures cannot be applied with success to 
every system of equations. However, when a convergent matrix 9TC 
can be found, the methods which have been discussed offer several 
distinct advantages over a straightforward matrix inversion. For ex- 
ample, if the order of the system is N, then it can be shown that each 
iteration requires approximately N 2 multiply-add operations.* On the 
other hand, an inversion requires at least N 3 multiply-adds, so that 
the relative saving is the ratio of the order N to the number of itera- 
tions required. It is often the case that the maximum eigenvalue is so 
small that the number of iterations required for an accuracy equivalent 
to that obtained by inversion is considerably less than N. 

Iterative methods also have the property that the solution accuracy 
is "adjustable", in the sense that once the solution has converged to 
the point where some norm, e.g., e E or e H defined previously, is less 
than a specified tolerance, the iteration process can be terminated. 
This property is especially attractive in view of the fact that truncation 
errors have already been introduced, and it would therefore be super- 
fluous to accurately invert a system which is itself approximate. By 
having the option of terminating the iterative procedure, we introduce 
an additional degree of freedom by which we can optimize the computa- 
tional program. 

V. APPLICATIONS 

The iterative techniques discussed in the previous section will now 
be applied to two problems of interest, namely the H-plane discontinuity 

* A multiply-add consists of the multiplication of two complex numbers and the 
adding of the residt to a third complex number. For repetitive computational 
algorithms, the number of multiply-ndds is a measure of the computational effort 
required. 
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problem mentioned in Section IV and the analysis of TE„ -* TM„ 
mode conversion at a step discontinuity in a circular waveguide. In 
both of these examples the required matrix elements may be expressed 
in convenient closed forms, which considerably reduces the required 
computational effort. 

5.1 H -plane Discontinuity in a Rectangular Guide 

In Section IV, the H-plane discontinuity of Fig. 3 was analyzed 
using an iterative algorithm which was observed to exhibit an oscillatory 
instability when initiated with a poor approximation to the actual 
solution. It was concluded that this was due to the eigenvalues of the 
iteration matrix 3TC being close to, or perhaps greater than unity. We 
now consider the same problem using the Gauss-Seidel method which, 
on the basis of the previous discussion, is expected to improve matters 
substantially. 

We assume a normalized TE 10 mode incident from the smaller guide. 
Because of the symmetry of the junction, such a wave excites only 
TE modes in both the + and — regions. The problem is, of course, 
to determine the corresponding modal coefficients for the fields on 
each side of the discontinuity. We shall present results only for the 
transmitted TE 10 mode, which is the only mode propagating in the 
larger waveguide for the present dimensions, ka t — 3.5, ka 2 = 4.5. 

It can be shown" that the normalized vector wave functions are 
given by 



'E'J = e. 



ico/x 



sin 



V2a^b -2a. 



pir 



(.1- + Oi] 



f E" = e t 



y _^sin[~g^(.x- + a 2 )~], 



(49) 



V2a,6 



where a, , a 2 , and b are the dimensions of the guide as shown and e„ is 
a unit vector in the y direction. The corresponding magnetic vector 
wave functions can be found from (12). The respective propagation 
constants are 

'*• - feT" 



~h" = 



*h" = 



2a J 



\i _ 

2 



(50) 



* - &T- 



From (21) and (49) we find that the coupling coefficients for the electric 
fields are given by 
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= 2, 



a { . pir . air 
zo) m \l~ sin ^r- sin ^r- 



?;tt + — gw 

<<2 



pr , d\ 7r 

""2 +^«2 



+-V- (>-£«!)}■ (51) 



«7T — — 07T 

a 2 



The appropriate magnetic field coupling coefficients are easily found 
from (51) using the relation 



3C, 



+ , - 



*Z. 



11 -rill~ 



", + 



where Z v denotes the modal impedance, equal to 



*Z" = 



con 

k h"' 



(52) 



(53) 



The results for the TE 10 modal coefficient, *B t , as obtained using 
the Gauss-Seidel iteration method, are conveniently represented in 
Table I. Also given are the numerical values for the error parameters 
e p , e B and £„ . Truncation for this example was at 25 modes in each 
waveguide. 

We conclude that for most applications, two iterations would probably 
have been sufficient, corresponding to a saving of greater than 90 percent 
in actual execution time, compared to a matrix inversion. The reason 
for this extremely rapid convergence is, as expected, in the magnitude 
of the largest eigenvalue, which was found, using (43), to be less than 
0.078. 

As a means of establishing the convergence of the normal mode 
solution, we have plotted in Figs. 7 and 8, the mean square errors 
£ E and i„ , respectively, as a function of the number of modes taken. 



Table 1^ — Results Using the Gauss-Seidel Method for Analysis 
of the H-Plane Discontinuity 



Iteration 
number 




li, 


Energy 
coefficient 


rma error 
<p 




Real 


Imaginary 


';/ 


1 
2 

3 
4 
5 


0.97445 
0.97747 
0.97747 
0.97747 
0.97747 


0.00951 
0.00455 
0.0042G 
0.00424 
0.00424 


-0.62X10- 2 

-0.74X10- 6 

0.85X10-' 

0.71X10- 6 

0.51X10-' 


0.5X10" 5 
0.49xnr B 
0.49X10" 5 
0.49X10- 5 
0.49X10-6 


0.01445 
0.01355 
0.01350 
0.01349 
0.01349 
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Fig. 7 — Mean square error e E , as function of the number of modes used, for 
H-plane discontinuity. 

These figures give genuine significance to the term "convergence in 
mean square", since they indicate that the use of more terms leads 
to better agreement with the boundary conditions in the mean square 
sense. Again, the error is uniformly higher for the magnetic field than 
for the electric field, due to the singularity at the corner of the dis- 
continuity. 

It is finally of interest to determine the effect of a poor initial estimate 
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Fig. 8 — Mean square error e H as function of the number of modes used, for 
//-plane discontinuity. 
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Table II — Results of Spoiling Electric Field at 4th Iteration 



Iteration 
number 


TEio coefficient 


Real 


Imaginary 


3 

4 
5 
6 

7 


0.97747 
1.94165 
0.98052 
0.97747 
0.97746 


0.00426 
0.02572 
-0.00063 
0.00396 
0.00423 



of the solution. We find that, unlike the simple algorithm discussed 
in the previous section, the Gauss-Seidel procedure is very stable, 
returning to the correct "steady state" solution within a few iterations 
after the spoiling was introduced. The results are shown in Table II, 
again for a truncation of 25 modes in each waveguide. 

5.2 Mode Conversion at a Step Discontinuity in a Circular Waveguide 

The second problem to which these techniques were applied is that 
of calculating the TE,, —* TM n mode conversion at an abrupt dis- 
continuity in a circular waveguide. Recent studies have indicated that 
this configuration is a very efficient transducer for use in dual mode 
conical horns. 10 The discontinuity is illustrated in Fig. 9 which shows 
the TE U mode, incident from the smaller guide, being converted to 
a combination of TE„ and TM U modes propagating in the larger guide. 

The normalized TE and TM vector wave functions are known" 
and, fortunately, it is possible to determine the appropriate coupling 
coefficients. For the elements of the matrix we find .that 




Fig. 9 — TEu — TMn mode conversion at a discontinuity in a circular waveguide. 
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b (a 



+ 



(K yl - yl) V(yl - i)(y 2 - l) /,(y.) 



*,/,(*.) Vfc£- 1) 



-, + 



= 0. 



(54) 



(55) 



(56) 



(57) 



In (59) through (62) we have used the following notation: 

a — radius of smaller waveguide, 

6 — radius of larger waveguide, 
x„ — pth zero of Ji(x), 
y p — pth zero of J[(y). 

The elements of the matrix 3C can easily be found from the impedance 
relations of (11) and (12). 

One parameter which has been found to be useful in characterizing 
the mode conversion properties of the discontinuity is the conversion 
coefficient C, defined as the ratio of the p-components of electric field 
for the two modes evaluated at the wall of the larger waveguide, i.e., 



C = 20 log,,, 



El 



E a 



dB. 



(58) 



p-6 



This quantity was calculated for the particular discontinuity a = 1.05", 
b = 1.4" over the frequency range 5.2 to 7.0 kHz, these parameters 
having been chosen for purposes of comparison with available experi- 
mental data. Truncation of the normal mode expansion was made 
after twenty-five TE and twenty-five TM modes in each waveguide. 
The iterative sequence was terminated when successive values of the 
modal amplitudes differed by less than 10" 6 . It was found that typical 
values for the error criteria are e P tt 10" r , and e E , £u ~ 0.015. These 
results indicate that for a given accuracy, a much lower value can be 
expected for e P , which is a function only of the lower-order modes, 
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Fig. 10 — TEn — > TMn conversion coefficient of a step discontinuity in a circular 
waveguide, a = 1.05", b = 1.4". 

than for e B and e H , which depend on the higher-order terms as well. 
In Fig. 10 we have plotted the computed values of the conversion 
coefficient, defined in (58), as a function of frequency. Also shown, 
us discrete points, are experimental results obtained previously. 10 The 
theoretical values are seen to be in very good agreement. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we have considered the solution of those matrix equa- 
tions which arise in the analysis of a class of waveguide discontinuity 
problems. In searching for criteria to estimate the accuracy of computed 
results, we have found that the uniqueness theorem itself yields a 
convenient set of error parameters which are easily implemented in 
the computational program. 

It is suggested that an iterative technique, particularly the "Gauss- 
Seidel" or "point-single-step" method often leads to a rapidly con- 
verging solution, thus offering several advantages over the usual inver- 
tive procedures, particularly in the speed of execution. Of particular 
interest is the fact that when this method is applied to the analysis 
of TE n — * TM„ mode conversion at an abrupt discontinuity in a 
circular waveguide, it yields a rapidly convergent and accurate solution. 
This has been established not only on the basis of theoretical error 
criteria, but also by comparison with experimental results previously 
obtained. 
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